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The van der Waals heterostructures of two-dimensional (2D) atomic crystals constitute a new paradigm in 
nanoscience. Hybrid devices of graphene with insulating 2D hexagonal boron nitride (h-BN) have emerged 
as promising nanoelectronic architectures through demonstrations of ultrahigh electron mobilities and 
charge-based tunnel transistors. Here, we expand the functional horizon of such 2D materials 
demonstrating the quantum tunneling of spin polarized electrons through atomic planes of CVD grown 
h-BN. We report excellent tunneling behavior of h-BN layers together with tunnel spin injection and 
transport in graphene using ferromagnet/h-BN contacts. Employing h-BN tunnel contacts, we observe 
enhancements in both spin signal amplitude and lifetime by an order of magnitude. We demonstrate spin 
transport and precession over micrometer- scale distances with spin lifetime up to 0.46 nanosecond. Our 
results and complementary magnetoresistance calculations illustrate that CVD h-BN tunnel barrier 
provides a reliable, reproducible and alternative approach to address the conductivity mismatch problem for 
spin injection into graphene. 

The discovery of graphene^'^ and several other 2D materials^ has opened up the possibility of employing their 
fascinating multifunctional properties for nanoelectronic and spintronic applications^"^. Graphene in par- 
ticular is an ideal material for spintronics^ due its predicted high electron spin coherence time and diffusion 
length^'^^. Since the first report on spin transport and precession in graphene^ \ vigorous studies are being pursued 
to achieve these theoretical goals^^'^l Electrical spin injection into graphene from a ferromagnetic (EM) metal is 
achieved by employing an insulating tunnel barrier, which circumvents the intrinsic conductivity mismatch 
problem^^. Inspired by the famous tunnel magnetoresistance structures, ferromagnetic tunnel contacts with 
metal-oxide barriers are often used for spin injection and detection in graphene^^'^^. The conventional metal- 
oxide barriers such as AI2O3 and MgO grown on graphene suffer from pin-holes, interface roughness and defects 
related to oxygen vacancies and doping in graphene^^"^^. Even though these barriers offer a certain level of 
interface resistance, there has been a lack of reproducible evidence for tunneling transport. Although studies 
involving Ti02 seeded MgO barriers^^ showed tunneling behavior, such hybrid metal-oxide barriers still suffer 
from roughness, non-crystallinity and reproducibility issues^°. These factors greatly influence spin polarization 
and spin lifetimes in graphene^^"^^. Therefore, the quest is still on for a better spin tunnel barrier, a prime building 
block in determining the future course of graphene spintronics^\ 

Here, we explore the role of two dimensional hexagonal boron nitride (h-BN) as a tunnel barrier for spin 
injection into graphene. In recent years, there has been a growing interest in van der Waals heterostructures of 2D 
crystals such as graphene with h-BN. In fact, van der Waals heterostructures can lead to the observation of many 
interesting phenomena such as tunnel transistors, metal-insulator transition, itinerant magnetism, Wigner 
crystallization^^"^^. Hexagonal boron nitride is an insulating isomorph of graphene with a bandgap of 6 eV^^ 
and a lattice mismatch of only 1.7%^^. Being a 2D crystal, it is ideally devoid of dangling bonds, atomically smooth, 
and possesses high mechanical strength and thermal conductivity, as well as even higher chemical stability than 
graphene. The h-BN has been successfully used as a substrate to improve mobility^^'^^ and to achieve long distance 
spin transport in graphene^^. Also, it has recently been employed as a tunnel barrier in graphene-based transis- 
^Qj.g22-24 'pj^g latter experiments revealed the possibility to precisely control the resistance by varying the number 
of atomic layers. The 2D structure of the h-BN makes it appealing for spin tunnel barriers because of the absence 
of surface states, which rules out any charge or spin traps at the interface^^. Such 2D layers can also act as diffusion 
barriers, which protect the underlying graphene from possible doping during the deposition of ferromagnetic 
contacts and subsequent processing. 
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Figure 1 | Hexagonal boron nitride (h-BN) tunnel barrier for graphene spintronics. (a), Scheme showing spin injection into graphene through 
h-BN tunnel barrier, (b), Schematic representation of the injection of spin polarized electrons from a ferromagnet into graphene through an atomically 
thin h-BN/ferromagnet tunnel structure, (c), Atomic Force Microscopy image of a CVD h-BN layer transferred onto a Si/Si02 substrate, (d), Colored 
optical microscope image of a multiterminal spin transport device showing a heterostructure of a graphene (Gr.) flake covered by a h-BN tunnel barrier 
and ferromagnetic cobalt electrodes patterned by electron-beam lithography. 



Recent theoretical studies by Karpan et al. showed that the tunnel 
contacts with tailored resistance area product (RA) and large spin 
polarization values can be achieved by engineering the number of h- 
BN layers at the junction^^. Studies also predict tunnel magneto re- 
sistance (TMR) exceeding 100% and perfect spin filtering with h-BN 
tunnel barriers^^'^^. In a recent experimental attempt, Yamaguchi 
et al.^^ used exfohated single layer h-BN for spin injection into gra- 
phene with spin signal of a few milliohms (with spin polarization < 
2%) and spin lifetime —50 ps. The low resistive h-BN barrier used in 
their experiments lead to no considerable improvements when com- 
pared to ferromagnet/graphene transparent contacts^'*. Also, the 
exfoliated h-BN used in these studies are smaller in size, which com- 
plicate the fabrication scheme and is unsuitable for wafer-scale pro- 
cessing. Therefore, in order to establish the significance of h-BN 
tunnel barriers for spin injection, it is necessary to demonstrate 
possibilities of achieving higher spin signals and spin lifetimes in 
graphene. At the same time it is also important to explore the viability 
of CVD grown h-BN as tunnel barrier for both charge and spin based 
devices, which has prospect for wafer scale fabrication. We address 
both these challenges in our work by employing CVD grown h-BN as 
a tunnel barrier for spin injection into graphene. 

In this letter, we show reliable and reproducible tunneling beha- 
viour of large area CVD grown h-BN and demonstrate efficient spin 
injection, transport and precession in graphene. Furthermore, we 
report an enhancement of the spin signal and spin lifetime by using 



resistive h-BN tunnel barriers, which circumvents the conductivity 
mismatch problem. Our observations supported by calculations elu- 
cidate that h-BN tunnel barrier resistances can be tailored to achieve 
desirable magnetoresistance values. While the present work focuses 
on the graphene spintronic devices, our results are generic for a large 
class of h-BN based charge and spin tunnel devices. 

Results 

We illustrate the underlying concept and spin injection device struc- 
ture in Fig. la and lb, which integrates the heterostructures of 2D 
materials (Graphene and h-BN) with ferromagnetic contacts. The 
spin polarized current is injected through a Co/h-BN ferromagnetic 
tunnel junction, creating a non-equilibrium spin accumulation and 
splitting of the chemical potential (A|i) in graphene (Fig. la). The 
resulting spin polarization is detected nonlocally by another h-BN/ 
Co tunnel contact placed at a distance of typical spin diffusion length 
(~2 jam). Using such van der Waals heterostructures, we dem- 
onstrate tunneling characteristics of large area CVD h-BN, tunnel 
spin injection signal at room temperature with a micrometer-length 
spin transport and Hanle spin precession in graphene. 

Tunneling characterization of hexagonal boron nitride. We 

fabricated graphene spin transport devices using h-BN as tunnel 
barrier and ferromagnetic Co electrodes for spin injection and 
detection. We prepared graphene flakes on Si/Si02 substrate by 
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Figure 2 | Electrical characterization of the Co/h-BN tunnel contacts, (a), Circuit configurations for 2-terminal (2T), 3-terminal (3T) and nonlocal (NL) 
measurement (b), 2-terminal and 3-terminal IV characteristics of a Co/h-BN/graphene contact with fits to Brinkman Dynes and Rowell (BDR) model. 
The inset shows the BDR fit to the dl/dV of 2-terminal measurements (c), Temperature-dependence of h-BN tunnel resistance measured in 3-terminal 
geometry at 170 mV (Top) and zero bias (ZBR) (Bottom), (d), Resistance area product (RA) (at zero bias) of the injection tunnel contacts obtained on 
different devices. 



exfoHation from HOPG using the scotch tape method\ Subsequently, 
we transferred the CVD grown layer of h-BN using a wet transfer 
process (elaborated in Methods and Supplementary information). 
We have been able to achieve a ripple free transfer of CVD h-BN 
over large areas on our Si/Si02 chips. The h-BN was then patterned to 
regions well covering the graphene flakes (see Supplementary 
information Fig. SI and Fig. S2). The position of the Raman shift 
for the CVD h-BN has been reported to be —1369 cm"^ by the 
supplier (graphene supermarket), which matches well with the 
thickness based studies on exfoUated h-BN^^ We also performed 
AFM (atomic force microscopy) characterization revealing a 
minimum thickness of 5 A for the h-BN layer as shown in Fig. Ic, 
which is the effective thickness observed for single layer h-BN^^'^^'^^. 
In our AFM based analysis, we also found the thickness of the h-BN 
to be about —5-10 A (Supplementary information Fig. S3) with few 
regions having even higher thickness. It has been reported that the 
resistance area product (RA) of h-BN scales exponentially with the 
number of layers^^. Thus, the RA scaUng analysis is an efficient 
method to resolve the thickness of h-BN, especially when the 
number of layers lies between 1-3. In our experiments, we used 
AFM and resistance scaUng analysis to ascertain the atomically thin 
nature of the CVD h-BN. Because of the variation in mean thickness, 
we also found the RA of the Co/h-BN contacts to vary between 
devices. Nevertheless, we found that the contacts in a particular 
device have reasonably similar RA indicating a uniform local 
thickness of the h-BN layer. 

We fabricated devices with and without h-BN tunnel barriers. In 
the following description, we refer h-BNO as the device without any 



tunnel barrier, and h-BNl, h-BN2, h-BN3 and h-BN4 as devices with 
h-BN tunnel barriers in the order of increasing RA. The details of the 
device fabrication are described in the supplementary information. A 
representative optical micrograph (color coded with enhanced con- 
trast) of a device with a graphene flake covered by a patterned atom- 
ically thin layer of CVD h-BN and several ferromagnetic electrodes 
fabricated by electron-beam lithography is shown in Fig. Id. 

We performed the electrical characterization of the devices in 
different measurement geometries (Fig. 2a) to ascertain the tunnel- 
ing behavior of the Co/h-BN/graphene contacts. The I-V character- 
istics were measured in two-terminal (2T) and three-terminal (3T) 
configurations. The channel characterization was done by regular 4- 
probe measurements and the spin signal was measured in nonlocal 
(NL) geometry. Fig. 2b shows the 2T and 3T tunneling I-V char- 
acteristic curves of an h-BN tunnel contact with RA —10 kQ[im^. 
The IV characteristics fit well (shown by red lines in Fig. 2b) to 
Brinkman-Dynes-Rowell (BDR) modeP^ for an asymmetric barrier. 
The fittings provide effective barrier heights of (p — 1.57 eV for 2T 
and 1.49 eV for 3T. These values match with the barrier heights for 
h-BN from previous reports^^'^^ and are close to typical barrier 
heights for oxide tunnel barriers used in magnetic tunnel junctions 
(MTJs). The weak temperature dependence of the tunnel resistance is 
a robust criterion for a tunnel junction without pin holes and trapped 
states^^. Such weak temperature dependence for the h-BN tunnel 
contact at zero bias and 170 mV is displayed in Fig. 2c. In addition, 
we also tested metal/h-BN/metal junctions where we obtained tun- 
neling IV characteristic curves fitting well to the BDR model apart 
from their weak temperature dependence. We find a variation in the 
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Figure 3 | Spin transport in graphene at room temperature, (a), Nonlocal spin valve signal (Rnl = Vnl/I) measured at room temperature for h-BNO 
(device with no tunnel barrier), h-BNl and h-BN2 (devices with different h-BN tunnel resistances) using an injection current of I = +5 \iA (Inset: 
Nonlocal measurement configuration). Sweep directions of the magnetic field are indicated by blue (positive to negative field sweep) and red (negative to 
positive field sweep) colors. The magnetic configurations of the injection and detection electrodes are illustrated for both sweep directions in the plot for 
h-BNl (b), Nonlocal spin valve voltage signal (AV = Vp - Vap) as a function of injection bias current for device h-BN2 at temperatures 100 K, 200 K and 
300 K. (c), Enhancement in spin signal (AR^l) with increase in tunnel contact resistance (RA) in different devices. 



RA product from device to device, which can be primarily attributed 
to the variation in mean thickness of CVD h-BN (5-10 A) as 
observed in AFM measurements. Fig. 2d shows the different RA 
obtained in our devices with h-BN tunnel barrier (Co/h-BN/ 
Graphene) and without h-BN (Co/Graphene) tunnel barrier. We 
note that tunneling characteristic curves with weak temperature 
dependence have also been seen in our other devices (h-BNl, h- 
BN2, h-BN3 and h-BN4, for example I-V curves for h-BN2 is shown 
in Supplementary information Fig. S4 and Fig. S5). These facts point 
out that our h-BN tunnel barrier is of good quality, and devoid of 
pinholes and interface states. The RA of the h-BN tunnel contacts in 
our devices is found to be between 2-13 kQ jim^, where we measured 
the spin signal. These values correspond to a thickness in the range of 
1-2 layers of h-BN in light of the recent thickness dependent resist- 
ance studies^^ (Supplementary information Fig. S6). The values of the 
RA for the Co/h-BN contacts are orders of magnitude larger when 
compared to transparent Co contacts in device h-BNO (—400 Q 
|im^). In spite of the challenges associated with the growth and trans- 
fer of CVD h-BN, we have been able to observe the tunneling beha- 
vior in most of the contacts over a large scale, over a wide range of 
temperature. The graphene channel resistivity in our devices was 
found to be in the range of 2-6 kQ with linear I-V behavior. We 
observed characteristic gate dependent Dirac curves for the graphene 
channels of our devices with a mobility of 2000-3000 cm^ V"^ s"\ 
The details of the graphene characterization are presented in 
Supplementary information (Fig. S7). 

Spin transport and precession in Graphene with h-BN/Co tunnel 
contacts. Having established the tunneling behavior of the Co/h-BN 



contacts, we focus on spin dependent transport measurements 
carried out on our graphene devices. The measurements were 
performed in a nonlocal (NL) spin valve geometry, where the 
charge current path is isolated from the spin diffusion signal 
measured by the voltage probes, as shown in the inset of Fig. 3a. In 
a typical device, we inject current between the +1 and -I contacts, 
and measure the nonlocal voltage signal between the +V and -V 
contacts. The spins injected through the Co/h-BN contacts 
accumulate in the graphene channel and diffuse laterally, thereby 
being detected at the nonlocal voltage probes. The nonlocal 
resistance is recorded while the in-plane magnetic field is swept 
from a negative to a positive value, followed by a reverse sweep. 
As shown in Fig. 3a, the spin signal is characterized by a change 
in the nonlocal resistance (ARnl^AV^l/I)- Rnl changes it's 
sign for parallel and anti-parallel configuration of the ferroma- 
gnetic electrodes in accordance with the changing spin density 
(majority to minority and vice-versa) at the detector electrode^\ 
Ferromagnetic electrodes having different shape anisotropies 
(as shown in Fig. Id) were employed as injectors and detectors 
to achieve different magnetization switching fields. The measure- 
ments were performed with a DC current in the range of ± 10 jiA. 

The room temperature nonlocal spin valve measurements for 
devices with a h-BN tunnel barrier (h-BNl and h-BN2) and without 
a tunnel barrier (h-BNO) are shown in Fig. 3a. The zero bias resist- 
ance values of the devices are presented in Fig. 2d. We observed 
nonlocal spin valve switching from parallel to antiparallel configura- 
tions of the injector and detector electrodes in our devices (with the 
magnetic configuration shown for h-BNl in Fig. 3a). The additional 
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Figure 4 | Spin precession measurement in graphene. (a), Nonlocal Hanle spin precession signal measured for devices with different tunnel barrier 
resistances (h-BNO, h-BN2, and h-BN3) as a function of the perpendicular magnetic field, keeping magnetization of the ferromagnetic electrodes in 
parallel in-plane configuration at 100 K, 200 K and 300 K. (b), Top: Bias dependence of the spin lifetime for device h-BN2 at different temperatures. 
Bottom: Bias dependence for h-BN3 at 100 K. (c), Dependence of spin lifetime on the contact resistance extracted from h-BNO, h-BN2 and h-BN3 at 
100 K. 



switching in device h-BN2 in the low field range is due to spin current 
reaching the outer reference Co electrodes. Notably, the magnitude 
of the spin signal is enhanced from h-BNO to h-BNl, h-BN2 and h- 
BN3 with increasing contact RA. In all the devices, we observed that 
the spin signal depends linearly on the injection current. We note 
that we did not measure a spin signal in our h-BN4 device with large 
tunnel contact resistance. A representative bias dependence of the 
spin signal for h-BN2 at different temperatures is shown in Fig. 3b. 
This confirms that the spin accumulation in the graphene channel 
created by the tunnel spin injection from the Co/h-BN contact is in 
the linear response regime. The enhancement in the nonlocal spin 
valve resistance AR^l with contact resistance RA for different devices 
is shown in Fig. 3c. We note that the spin signal shows a monotonic 
increase in value as the h-BN tunnel contact resistance is increased. 
The insertion of a high resistive h-BN tunnel barrier not only pre- 
vents the backflow of spins into the ferromagnetic contact, but also 
minimizes the spin flip scattering at the interface leading to a higher 
spin accumulation in graphene and hence a higher spin signal. We 
describe later (in the nonlocal spin signal calculation) how the AR^l 
is quantitatively related to the tunnel barrier resistance. 

In order to evaluate the spin lifetime of the electrons in graphene 
devices with h-BN tunnel barrier, we performed nonlocal Hanle spin 
precession measurements by sweeping the magnetic field perpendic- 
ular to the device geometry while the magnetization axes of the 
injector and detector electrodes are being kept parallel. The injected 
spin polarized electrons precess about the perpendicular magnetic 
field with Larmor frequency (wl) while diffusing towards the non- 
local detector contact. In Fig. 4a, we show the nonlocal Hanle signal 
measured for devices with h-BN tunnel barrier (h-BN2 and h-BN3) 
and without tunnel barrier (h-BNO) at different temperatures for a 
bias current of + 5 jiA. In addition to the increase in the magnitude 
of the spin signal we also observe a decrease in line width of the Hanle 
curves. For device h-BN3 with contact RA = 12.6 |am^, we 
observed a nonlocal Hanle resistance with AR^l ~ 8.5 Q at 100 K. 
This would correspond to a nonlocal spin valve signal of 17 Q. To 
remove possible non-spin related effects, we have symmetrized our 



raw Hanle data by taking the average of the negative and the positive 
magnetic field sweeps^^. To estimate the spin lifetime, we fit the 
measured data using the Hanle spin transport equation (Eq. 1), 
involving the spin diffusion, precession, and relaxation factors. The 
variation of nonlocal resistance (ARnl) due to Larmor precession of 
the spins diffusing towards the detector from the injector is given by 

AR^L = + — ^"^ cos(coit) dt (1) 

- Wash ^4nDt ^ ' ^ ^ ^ 

where P is the spin polarization of the h-BN/Co contacts (assumed to 

be the same for both the injector and the detector), and col = B± 

n 

(with Lande's g-factor g = 2) is the Larmor frequency of the preces- 
sion of electrons in a perpendicular magnetic field Bj^. L is the 
effective channel length between the injector and the detector 
(2 jim in our devices), with a width W and square conductivity 
of the graphene strip. The diffusing spin polarized electrons have a 
lifetime Ts and diffusion constant D with a spin diffusion length 

The spin lifetime in graphene devices with h-BN tunnel barrier (h- 
BN2, h-BN3) are found to be much higher than the device without 
tunnel barrier (h-BNO). When compared to h-BNO having ~ 
50 ps, the evaluated are higher for devices h-BN2 (t^ ~ 110 ps 
with RA ~ 5.8 kQ jim") and h-BN3 (t, ^ 0.46 ns with RA = 
12.6 kQ |im^). The maximum spin signal and spin lifetime were 
obtained for the h-BN3 device, giving rise to a spin diffusion length 
Xsf of 1.4 jam with a spin polarization P of 14% and —0.46 ns. The 
diffusion constant obtained from the Hanle fitting is D = 
0.00416 m^/s. The spin lifetime T5 is found to be independent of 
the bias current for our devices as shown in Fig. 4b. We further 
observe that for h-BN2, the spin lifetime is also fairly independent 
of the temperature (top panel of Fig. 4b). In this case the bias current 
does not influence the diffusive spin transport in graphene much. 
The weak temperature dependence of the spin lifetime also indicates 
that the spin relaxation in the channel is mainly dominated by the 
impurity and charged scattering centers in graphene. Most impor- 
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Figure 5 | Calculated spin signal as a function of tunnel contact resistance, (a), Nonlocal resistance vs. tunnel barrier resistance with the experimental 
value marked by the red dot. Inset: nonlocal configuration, (b), Dependence of 2-terminal magnetoresistance (MR) on the tunnel barrier resistance with 
the arrow in the plot indicating the tunnel resistance (Rij) in device h-BN3. Inset: 2 terminal configuration. 



tantly, the magnitude of the spin signal and the spin lifetime are 
found to increase with tunnel barrier resistance. A plot showing 
the enhancement in the spin lifetime with barrier resistance is dis- 
played in Fig. 4c. Such behavior can be understood considering the 
various sources of scattering. The spin relaxation time (Xs) is prim- 
arily a measure of spin flip scattering (Xsf"^) inside the graphene 
channel. In the case of spin injection from a ferromagnet to graphene, 
the spin escape rate (Xesc^)j ^ rate by which the spins diffuse back 
into the ferromagnet and get relaxed, greatly influences the spin 
lifetime (t ~ ^ = t^^ ^ + t~J). A tunnel barrier overcomes this problem 
by reducing the backflow probability (ideally, Xesc -^^). In addition 
to this, the spin accumulation at the interface suffers from further 
relaxation due to spin flip scattering at the ferromagnet- graphene 
interface^^. The resistance of the tunnel barrier is a measure of the 
effective isolation to avoid backflow and spin flip scattering at the 
interface. Thus an increase in h-BN tunnel barrier resistance results 
in the reduction of contact induced spin relaxation^^'^^, leading to an 
enhancement in the spin signal and lifetime. However, an indefinite 
increase of barrier resistance is not viable and we discuss about the 
upper range of tunnel barrier resistance for maximum spin accu- 
mulation in the following. 

Nonlocal and local magnetoresistance calculations. To understand 
the effect of barrier resistance on the overall spin signal, we calculate 
the nonlocal spin signal as a function of the tunnel contact resistance 
(Rt) following the drift- diffusion based analytical formulation by 
Takahashi and Maekawa^°. Fig. 5a displays the variation of the 
spin signal with tunnel resistance for a graphene nonlocal spin 
valve, calculated using the obtained polarization and spin lifetime 
of h-BN3. The curve indicates that the nonlocal spin valve signal 
increases with increasing tunnel barrier resistance and saturates 
beyond an optimum value of the tunnel resistance (Rid = yjRiRd, Ri 
~ R^y /^injector, (i-^ detector). The value of the nonlocal spin valve 
signal ARnl for h-BN3 (~ twice the Hanle signal in h-BN3) is 
marked by a red dot along with the theoretical curve. The close 
agreement between the calculated spin signal (which is dependent 
on tunnel resistance) and the observed spin valve signal implies that 
the values of obtained polarization and spin lifetime are reliable. The 
plot also reveals that the tunnel resistance Ri^ in our device h-BN3 lies 
close to an optimum value for the observation of the maximum spin 
signal. Beyond the contact induced relaxation regime, any further 
increase in tunnel barrier resistance will not only saturate the spin 
signal but will also result in higher power dissipation in the device. In 



addition to the saturation of the spin signal, an enhancement in the 
tunnel barrier resistance inevitably increases the dwell time (t„) of 
spins inside the channel. This leads to a down turn (as t„ > T5) in 
magnetoresistance (MR) as observed in spin valves. In Fig. 5b, we 
show the 2-terminal MR calculation for a channel with the spin 
parameters of device h-BN3 using Fert-Jaffre formula^^. This 
calculation shows the possibility of observation of a significant MR 
in a narrow range of contact resistances with R^^ of h-BN3 being 
closer to maximum. In the calculation, the low MR for small 
contact resistances is due to the conductivity mismatch, preventing 
spin injection from the Co contact into graphene. The decay of the 
MR for large contact resistances is due to relaxation of the spins 
during the time spent in the graphene channel (dwell time), a time 
that becomes very long for large contact resistances. Note that the 
values of spin parameters such as the spin polarization might also 
depend on the tunnel barrier resistance {Rt) which was not taken 
into account in these theoretical models. In spite of the fact that this 
might lead to slight deviation on either side of Rr ~ Rid^ the constant 
values that we employed here are quite reliable in the regime where 
~ ^id evident from the match in the non-local calculation. This 
analysis elucidates that CVD grown h-BN is a reliable solution to the 
conductivity mismatch problem without requiring the deposition of 
any oxide based tunneling barriers on graphene. It should be noted 
that the optimum tunnel resistance (i^^^^) for the observation of a 
maximum spin signal/MR also depends on the spin resistance {Rch 
oc ksf) and thus the quality of the channel. This implies that channels 
supporting larger spin diffusion length or a higher spin lifetime can 
overcome the dwell time limit set by the tunnel barrier resistance. 
With epitaxial graphene channels, Dlubak et al.^^ have shown large 
spin signal values with a high contact resistance, which has been 
attributed to the better quality of the epitaxial graphene. Therefore, 
it might be possible to observe higher MR values with better quality 
graphene while using h-BN layers as tunnel barriers having higher 
RA values. Furthermore, with the improvements in the h-BN CVD 
growth process, it should be possible to achieve precise control over 
the RA of the tunnel contacts. 

Discussion 

The current study using CVD h-BN tunnel barrier has several impli- 
cations. CVD grown 2D crystals hold the key for wafer scale proces- 
sing and technological implementation. Often, it has been a 
challenge to observe performance in CVD grown materials even 
comparable to their exfoliated counterparts, mostly due to the imper- 



SCIENTIFIC REPORT: | 4 : 6146 | DOI: 1 0.1 038/srep06146 



6 



fection introduced during their growth and subsequent transfer pro- 
cesses. In spite of these challenges, we have been able to observe 
tunneling behavior of CVD h-BN in all devices with h-BN coverage 
over a large area. On the other hand, most of the metal oxide barriers 
used in this field do not provide reproducible tunneling nature except 
for a few reports^^'^\ The tunneling characteristics are reproducible 
in our CVD h-BN devices despite the observed variation in the 
contact RA due to the spatial variation in thickness, which we have 
exploited to investigate the enhancement of spin signals as a function 
of the tunnel resistance. 

The values of spin lifetime of 0.46 ns and polarization of 14% 
obtained in our devices show that h-BN is a promising tunnel barrier 
for spin injection into graphene. This performance is more than an 
order of magnitude higher than the recent report by Yamaguchi et al. 
using exfoliated monolayer h-BN as tunnel barrier. The spin para- 
meters using CVD h-BN are also comparable to that observed using 
conventional MgO and AI2O3 barriers^ The widely observed spin 
lifetimes vary between 100-500 ps^^, except for occasional reports of 
spin lifetimes up to 1 ns and above in pristine and hydrogenated 
graphene^^'^^'^^ The fact that CVD h-BN can serve as a tunnel barrier 
capable of overcoming the conductivity mismatch problem over a 
large scale together with its potential makes it highly favorable. Due 
to the 2D nature of h-BN, it is highly compatible with graphene and 
bypasses the disadvantages of conventional metal- oxide barriers 
related to growth and doping in graphene. 

In addition to the quality of the tunnel barrier, the spin lifetime 
and polarization depend to a great extent also on the intrinsic scat- 
tering mechanisms inside the graphene medium. The observed spin 
lifetime and diffusion length are significant for our graphene flakes 
on Si02 substrate with a carrier mobility of 2000-3000 cm^ V"^ s~\ 
which can be enhanced further by improving the mobility of the 
graphene^ This, combined with a better quality of CVD h-BN 
certainly holds the key to answer the important questions in gra- 
phene spintronics. The prospect of atomically thin h-BN for tunnel 
spin injection goes well beyond this demonstration. Smooth and 
atomically flat h-BN tunnel barriers are expected to show higher spin 
lifetimes and spin injection efficiencies in graphene than their metal- 
oxide counterparts. This is highly intriguing, and encourages more 
detailed spin transport studies with h-BN tunnel barriers on high 
mobility graphene devices such as h-BN encapsulated graphene or 
suspended graphene channels with boron nitride support under the 
contacts. Recent studies elucidated the use of 2D graphene barrier for 
spin injection in silicon'*^ and magnetic tunnel junctions^^'^^. 
However, graphene, being a semi-metal with zero bandgap, is a 
low resistive barrier. In contrast, h-BN is a high bandgap insulator 
with a similar structure and can support spin injection into materials 
with larger spin resistances (Rch) such as graphene itself Oppor- 
tunities are growing with the proposal by Karpan et al.^^, showing 
layers of h-BN can be combined to the perfect spin filtering prop- 
erties of Co/MLG (multilayer graphene) to realize Co/MLG/h-BN 
ideal spin injectors. In addition, tunnel contacts with tailored RA and 
spin polarization values can be achieved by engineering the number 
of graphene and h-BN layers at the junction^^. Such contact RA 
engineering with improved spin parameters can potentially enhance 
magnetoresistance in graphene tunnel junctions^^'^^. Such hetero- 
structures of graphene and h-BN have a great potential for spintro- 
nics in general and are certainly systems for more interesting studies 
in the future. 

In conclusion, we have demonstrated tunnel spin injection into 
graphene using atomically thin CVD hexagonal boron nitride (h- 
BN) at room temperature. The CVD h-BN showed reproducible 
tunneling characteristics with weak temperature dependence of the 
tunnel resistance over a large scale. Taking the advantage of the 
variation in the thickness of CVD h-BN, we have been able to invest- 
igate the evolution of the spin transport with resistance of the h-BN 
tunnel contacts. With resistive h-BN tunnel contacts, we observed an 



enhancement in spin transport with spin lifetimes up to 0.46 ns, 
yielding a spin polarization of 14%. In light of the magnetoresistance 
calculations, we show that the CVD h-BN is a reliable and unique 
way of creating tunnel resistances, which overcomes the conductivity 
mismatch problem. The integration of CVD h-BN to graphene opens 
up new avenues to explore the role of h-BN atomic planes for spin 
injection into other 2D advanced materials^'^^, metals^^, semi- 
conductors^°'^\ and magnetic tunnel j unctions possibly leading 
to the observation of large magnetoresistance effects. 

Methods 

Graphene preparation. The graphene flakes were prepared on a clean Si02 
(285 nm), highly doped silicon substrate with predefined Ti/Au markers. We used 
graphene exfoliated from highly oriented pyrolytic graphite (obtained from 
Advanced Ceramics) using the conventional scotch tape method. Graphene flakes 
were identified using a combination of optical and atomic- force microscopy. In our 
experiments, we used graphene flakes with widths around 1.5 [im. 

h-BN tunnel barrier. The h-BN monolayers used as tunnel barriers were grown on 
copper substrate and have high crystalline quality with ~ 95% coverage (obtained 
from Graphene supermarket). During CVD growth, h-BN grows on both sides of the 
copper foil (h-BN/Cu/h-BN). The h-BN on the top surface was first protected with 
thick PMMA (Al 1 e-beam resist), and the h-BN on the lower surface was removed by 
argon plasma etching. Then the PMMA/h-BN/Cu structure was kept in FeCls 
solution for etching the Cu substrate. After complete removal of Cu, the floating layer 
of PMMA/h-BN film was washed several times with deionized (DI) water. To remove 
possible chemical residues, the film was also treated with 10% HCl for 20 min, 
followed by neutralization with DI water. The PMMA/h-BN film was then 
transferred onto graphene flakes prepared on a Si02/Si chip. The chip containing the 
PMMA/h-BN/graphene heterostructure was dried and baked at 150°C for 5 minutes 
to improve the adhesion between the h-BN and the graphene. Finally, the PMMA 
layer was removed using hot acetone and the h-BN/graphene on the chip was 
annealed at 200°C in an H2/Ar atmosphere for 20 minutes. Careful transfer processes 
with thicker PMMA and the subsequent optimized baking procedure ensured a 
uniform and ripple-free h-BN layer on the graphene. The h-BN layer was patterned to 
shapes covering the graphene flakes using electron-beam lithography and argon 
plasma etching (shown in Fig. 2a). 

Ferromagnetic contacts. Ferromagnetic electrodes were fabricated on h-BN/ 
graphene samples by electron-beam lithography and lift-off techniques. Cobalt thin 
film contacts of 65 nm thickness with a cap-layer of 15 nm gold were prepared using 
electron-beam evaporation. Ferromagnetic electrodes of different widths (130 nm, 
260 nm, 500 nm; aspect ratio > 15) for the spin injector and detector were fabricated, 
ensuring different switching fields to achieve a parallel or antiparallel configuration 
through in-plane magnetic field sweeps. Most of the data presented here correspond 
to contacts with 130 nm, 250 nm as the injector and detector electrodes. 

Electrical measurements. Electrical characterization of the contacts and the device 
were performed using multiterminal measurement geometry. The contact resistances 
were found to be in the range of 10-100 kQ, showing reasonable scaling with area, 
thus demonstrating the homogeneity of the h-BN tunnel barrier. The high contact 
resistance and its low temperature dependence (shown in Fig. 2c) confirm the high 
quality of the h-BN tunnel barriers on graphene. The channel resistance of the 
graphene measured using the four-probe technique was found to be 2-6 kQ. The 
characterization of graphene with gate shows Dirac curves with Vd~ 10-20 V and 
mobility ~ 2000-3000 cm^ V"^ s~\ 

Spin transport measurements. Spin transport and spin precession measurements 
were performed in nonlocal geometry where the charge current path is separated 
from the voltage contacts measuring spin signal, in order to exclude spurious signals 
such as anisotropic magnetoresistance and stray Hall voltage. The measurements 
were performed using direct currents (DC) of 1-5 |iA, and the nonlocal voltage was 
measured using a nanovoltmeter, under vacuum in a variable temperature cryostat 
(from room temperature 300 K to 1.5 K) with a superconducting magnet. 
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